Measuring the one-particle excitations of ultracold fermionic atoms 
by stimulated Raman spectroscopy 
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We propose a Raman spectroscopy technique which is able to probe the one-particle Green's 
function, the Fermi surface, and the quasiparticles of a gas of strongly interacting ultracold atoms. 
We give quantitative examples of experimentally accessible spectra. The efficiency of the method 
is validated by means of simulated images for the case of a usual Fermi liquid as well as for more 
exotic states: specific signatures of e.g. a d-wave pseudo-gap are clearly visible. 
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The remarkable advances in handling ultra-cold atomic 
gases have given birth to the new field of "condensed 
matter physics with light and atoms" . Key issues in the 
physics of strongly correlated quantum systems can be 
addressed from a new perspective in this context. The 
observation of the Mott transition of bosons in optical 
lattices pQ, of the superfluidity of fermionic gases [2], and 
the recent imaging of Fermi surfaces [3] have been im- 
portant milestones in this respect. Ultimately fermionic 
atoms in optical lattices [U |S] could help understanding 
some outstanding problems of condensed matter physics, 
such as high-temperature superconductivity. In this con- 
text, a key issue is the nature of the low-energy ex- 
citations of low-dimensional strongly interacting Fermi 
systems. There is abundant experimental evidence that 
those are highly unconventional, departing from standard 
Fermi-liquid theory. 

In this letter, we study how to probe the one-particle 
excitations of interacting ultracold fermionic atoms us- 
ing stimulated Raman spectroscopy. This technique 
has been considered previously in the context of cold 
atomic gases, as an outcoupling technique to produce 
an atom laser [BJ, and also as a measurement technique 
for bosons EJ QD] and fermions pU E2]- Here, we 
demonstrate that this technique provides, for strongly- 
interacting fermion g 3SGS ■ cL momentum-resolved access 
to key properties of the quasiparticle excitations, such as 
their dispersion relation and lifetime. It also allows for a 
determination of the Fermi surface itself in strongly inter- 
acting regimes, whereas previously demonstrated meth- 
ods [3] apply to the non-interacting case. Furthermore, 
it is shown that the suppression of quasiparticles due to 
a pseudogap in the excitation spectrum can also be de- 
tected by this method. 

In a conventional Fermi liquid, low-energy excitations 
are built out of quasiparticles [13] . Those are character- 
ized by their dispersion relation, i.e the energy £k (mea- 
sured from the ground-state energy) necessary to create 



such an excitation with (quasi-) momentum k. The in- 
teracting system possesses a Fermi surface (FS) defined 
by the location in momentum space on which the ex- 
citation energy vanishes: £k F — 0. Close to a given 
point on the FS, the quasiparticle energy vanishes as: 
£k ~ v>(kp) • (k — kp) + • • • , with vp the local Fermi 
velocity (inversely related to the effective mass). Quasi- 
particle excitations have a finite lifetime r^ 1 , and are well 
defined provided vanishes faster than £k as the FS is 
approached (Fk ~ in Fermi liquid theory). In con- 
trast, one-particle excitations in the "normal" (i.e non- 
superconducting) state of the cuprate superconductors 
(SC) reveal strong deviations from this behaviour [14] . 
Reasonably well-defined quasiparticle excitations only 
exist close to the diagonal direction of the Brillouin zone 
(the "nodal" direction along which the e?-wave gap van- 
ishes in the SC phase), and even there Tk is rather large. 
Away from this direction (in the "antinodal" region) , ex- 
citations appear to be short-lived and gapped already 
above the SC critical temperature (the so-called pseudo- 
gap phenomenon) . This momentum-space differentiation 
is a key to the physics of cuprates. 

Experiments probing directly non-diagonal one- 
particle correlators (ip^(r, t)ijj{r' , t')) of a many-body 
system arc therefore highly desirable but also rela- 
tively scarce. Most physical measurements indeed pro- 
vide information on two-particle correlators of the form 
(ipi(r,t)i/)(r,t)^(r',t')ip(r',t')) [24. Examples are neu- 
tron scattering or transport measurements in the solid- 
state context [13] , and Bragg scattering [15] or noise cor- 
relations measurements [16 in the cold atom context. 
For Bose systems with a finite condensate density uq, 
the two-particle correlator is closely related to the one- 
particle correlator via terms such as tlq (tp'(r, t) ip(r', t')). 
By contrast, in Fermi systems, the distinction between 
one- and two-particle correlators is essential and specific 
measurement techniques of the former are requested. 

In solids, angle-resolved photoemission spectroscopy 
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(ARPES) provides a direct probe of the one-particle 
spectrum |17j . and has played a key role in revealing 
momentum-space differentiation in cuprates [14]. It con- 
sists in measuring the energy and momentum of electrons 
emitted out of the solid exposed to an incident photon 
beam. In the simplest approximation, the emitted inten- 
sity can be related to the single-electron spectral func- 
tion, defined at T = and for ui < 0, i.e for hole-like ex- 
citations by: A(k, u) = £„ \{^- x \^\^)\ 2 6(u + M + 
E n — E ). In this expression, ip k is a destruction operator 
for an electron with momentum k, W q is the ground-state 
of the N-particle system and v I'^ r_1 are the eigenstates of 
the system with one less particle. In a conventional Fermi 
liquid, and for momenta close to the FS, the spectral 
function can be separated |13j into a coherent quasipar- 
ticle contribution and an incoherent contribution : A = 

A QP + A iac , with 7iviQp(k, uj) ~ z k r k /[(w - a) 2 + r k ] 

and Ai nc widely spread in frequency. Only a finite frac- 
tion Z k < 1 of the total spectral weight corresponds to 
long-lived quasiparticle excitations. 

In this paper, we consider stimulated Raman spec- 
troscopy on a two-component mixture of ultracold 
fcrmionic atoms in two internal states a and a' . Atoms 
are transferred from a into another internal state (3 ^ 
a, a', through an intermediate excited state 7, using two 
laser beams of wavectors ki : 2 and frequencies wi i2 . If 
lj\ is sufficiently far from single photon resonance to the 
excited 7 state, we can neglect spontaneous emission. 
Eliminating the excited state, we write an effective hamil- 
tonian: V = C J dr ^(r)V>/3(r)e i(kl_k2) ' r a|a 2 + h.c, in 
which a\ (02) denotes the creation (destruction) operator 
of a photon respectively in mode 1 (2) and the constant 
C is proportional to the product of the dipole matrix ele- 
ments d a ~[ and dpy of the optical transitions and inversely 
proportional to the detuning from the excited state. 

The total transfer rate to state (3 can be calculated [7J 
|8j[9] using the Fermi golden rule: 

i2(q,n) = |a| a ni(na + l) f dt f dr dr' e *[nt-q-(r-r')] 

J — OO J 

x 5/3 (r,r';i)(^+(r,t)^(r',0)). (1) 

Here q = ki — k 2 and O — u\ — L02 + ^ with /i the chemical 
potential of the interacting gas, and ni,2 the photon num- 
bers present in the laser beams. Assuming that no atoms 
are initially present in (3 and that the scattered atoms 
in /3 do not interact with the atoms in the initial a, a' 
states, the free propagator for /3-state atoms in vacuum is 
to be taken: g (r,r';t) = (0/jI V/j(r,t)Vj(r / s 0)|0 j9 ). The 
correlation function entering (JTJ) is proportional to the 
one-particle Green's function [25] (ip^(r,t)tp a (r',0)) = 
—iGa(r',r,—t) of the strongly- interacting Fermi sys- 
tem. For a uniform system, the rate (JTJ can be related 
to the spectral function A(k, ui) of atoms in the internal 



state a by [S]: 

i?(q,0)cx Jdkn F (e kp -n)A(k-q,£ kl3 -n) (2) 

in which the Green's function has been expressed in terms 
of the spectral function and the Fermi factor np as [13] : 
G< (k, to) = i n F (uj) A(k, u). 

In order to physically understand which information 
can be extracted from a measurement of the rate ^ , let 
us first approximate the spectral function by A(k,oj) = 
S(lu — £ k ), i.e neglect the incoherent part and consider 
quasiparticles with an infinite lifetime. The Raman rate 
then reads at T = 0: R = J^ k<Q dk <5(£k+ q ,/3 — £k — 
Contributions to this integral come from momenta in- 
side the FS (£t < 0) which satisfy the Raman reso- 
nance condition £k+q,/3 — £k = ^- When the frequency 
shift is small, R vanishes since there is no available 
phase-space satisfying these constraints. The smallest 
frequency at which a signal starts to be observed is 
Qt = Mink£k/3 = ejg [2B]. This corresponds to a mo- 
mentum transfer q = — k^ which lies itself on the FS 
(i-e £k F = 0) [27 ]. For very close to the extinction 
threshold (AO = il — £It ^ 0), the region in momentum 
space inside which a sizeable transfer rate R is measured 
consists of a shell surrounding the FS, centered at a mo- 
mentum q such that AO = — £_ q ~ vpfap) ■ (q + k^), 
and of width Aqu ~ v2M AO. In these expressions, M 
is the effective mass at the bottom of the /3-band and Vp 
is the Fermi velocity. 

This analysis remains unchanged when considering 
quasiparticles with a finite lifetime (uniform along 
the FS), the width of the momentum shell being simply 
replaced by Aqu ~ y2M AO + T/vp. Hence, measuring 
the Raman signal for O close to the extinction thresh- 
old Ot and sweeping over q, provides a determination of 
the FS in an interacting system (while the method of [3] 
applies to non-interacting fermions). It also gives access 
to the velocity of quasiparticles (from the displacement 
of the measured signal as a function of O) and to their 
lifetime (from the width of the momentum-shell) . 

Examples of numerically simulated spectra for uniform 
interacting systems are given in Fig. [lji-b, where a color 
intensity plot of the Raman rate ^ is shown for a fixed 
value of the frequency shift close to threshold. In (a), 
we consider the case of a Lorentzian spectral function 
centered around the free dispersion relation of a two- 
dimensional square lattice: £k = — 2< Q (cos fc^+cos k y )— /x. 
In (b) a phenomenological form [18 of the spectral func- 
tion is used, which captures the main aspects of the 
ARPES data in the non-SC ("normal") state of high- 
temperature superconductors. The key feature entering 
this phenomenological form is a pseudo-gap with <i-wave 
symmetry Ak = Ao(cosfc K — cosk y ), corresponding to a 
depletion of low-energy excitations even when no long- 
range SC order is present. Ak vanishes along the zone 
diagonal (nodes) and is maximum along (0, 0) — (w, 0) 
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FIG. 1: Intensity plots of the Raman rate -R(q, SI), for SI close 
to threshold Sir (ASl = 0.01t a ). (a) Non-interacting fermions 
on the homogeneous 2D square lattice with density n a = 0.22 
and a Lorentzian broadening of the spectral function V = 
QAt a uniform in k-space. (b) Model d-wave pseudogap state 
(see text), with A = 0.1 t a , T = 0.05t Q ,ri = 0At a . The 
plot is for a hole-doped system (n a — 0.45) with a nearest (t a ) 
and next-nearest neighbour (t' a ) hopping, with t' a /t a = —0.3 
(typical for cuprates, but similar effects are expected also for 
smaller \t' a /t a \). (c) and (d): same as (a) and (b) in the 
presence of a harmonic trap (uo = 0.02t a ). The pseudogap 
and nodal-antinodal differentiation are clearly visible in both 
lb and Id. 

(antinodes). A self-energy A^/(w + £k) is a conve- 
nient modclization of this effect, fn addition, finite 
lifetimes effects are introduced, resulting in the form: 
7rA(k, u) = -1m [« - £ k + iTi - A£ /(« + £ k + iT )] ~\ 
This corresponds to a quasiparticle dispersion which is 
gapped out except at the nodes. The width Tk = 
ri + A£r /[(w + £k) 2 +rg] is largest near the antinodes. 
This form of ^4(k, oS) also provides a reasonable qualita- 
tive description of recent theoretical results for the two- 
dimensional Hubbard model |19j . The momentum space 
differentiation encoded in the model spectral function is 
clearly visible on Fig. [TJd, with nodal regions display- 
ing quasiparticles while antinodal ones are gapped out 
and short-lived. This illustrates how the Raman spec- 
troscopy method can be used to determine the FS not 
only of a Fermi liquid but also of a strongly interacting 
system with suppressed quasiparticles. In Fig. [2] we fur- 
ther show that the lineshape of the Raman signal for a 
fixed value of q does reveal the essential features of the 
spectral function, namely quasiparticles at the nodes and 
a pseudo-gap at the antinodes. 

Since most cold atom experiments are performed in 
a trap, it is important to verify that the spatial in- 
homogeneity does not spoil the predicted signatures. 
Within the local density approximation, the observed 
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FIG. 2: Comparison between the spectral function A and the 
Raman rate R/ASl for two points in momentum-space indi- 
cated on Fig. 1. In the nodal direction (N; left), the spectrum 
displays a quasiparticle peak, while in the antinodal direction 
(A; right) a depletion of the signal is observed at low energy, 
corresponding to the pseudogap. 



signal is the sum of the contributions of the different 
points R of the trap, with a local chemical potential 
= A* — A/(x>qR 2 /2. The results are summarized in 
Figs. [ljc,ci) for physical situations such that the value of 
the chemical potential at the trap center coincides with 
that of the homogeneous system in Figsfija,b). As ex- 
pected, the intensity map is now a superposition of the 
Fermi surfaces corresponding to all the densities realized 
in the trap. The outer shell delimited by the extinction 
of the signal still gives a direct access to the FS cor- 
responding to the highest densities at the centre of the 
trap. The typical signatures of an unconventional state 
remain clearly visible in the trap as well: in Fig. [TJi, the 
nodal-antinodal differentiation is apparent in the outer 
shell of this plot, as seen from the suppressed intensity 
along the antinodal direction. A possible way of reveal- 
ing the region around the Fermi surface is to measure 
the intensity maps for two, slightly different values of the 
frequency and/or the total atom number, and then take 
their difference: the resulting differential images for the 
trapped system (not shown) recover the same qualitative 
features of the homogeneous system shown in Fig. [T^L-b. 

The discussion so far has assumed that it is possible 
to repeat the measurement of the total rate R for several 
different values of q. In some cases, a different scheme 
with a momentum-selective detection of the scattered /3 
atoms may be instead favorable, quite similar to ARPES 
in solids. A single value of q is used, and a time of flight 
expansion of the f3 atoms cloud is performed (after sud- 
denly turning off the trap and the lattice potential) in 
order to reconstruct the momentum distribution of the 
atoms. As shown in Fig. [3^, the Raman resonance con- 
dition allows for a selective addressing of the different 
regions in k by tuning the frequency SI. The number of 
Raman-scattered atoms with final momentum k is pro- 
portional to the integrand np^s^p — SI) A(k — q, £k/3 — ^) 
of Fig. 3b shows that the resulting k-space intensity 
map is able to reveal the details of the pseudogap physics, 
in particular its k-dependence. By varying both SI and 
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FIG. 3: (a) Colour plot illustrating the selective addressing of 
k-space by a proper choice of (cf. colour scale on the left), 
(b) Time of flight k-map obtained by integrating the Raman 
intensity for Afl = fi — fir varied in the range [2At a , 6.8t a ]. 
The dispersion relation of the /3-atoms is taken as £k/3 = £/j — 
2tp(2 + cosk x +cosk y ) with tp = 1.5t Q (note that interactions 
will renormalize downwards the effective t a even if bare values 
are equal). Parameters are as in Fig.lc and q = 0. 



q, Raman scattering offers more possibilities for probing 
the system in a momentum-selective way than microwave 
spectroscopy techniques [2D] . 

As a final point, we discuss some orders of magni- 
tude which are important for the actual feasibility of 
the experiments proposed in this article. Specifically, 
we consider 6 Li atoms (see also Ref. [12]) in the two 
lowest hyperfine states \a) = \I Z = l,m = —1/2) and 
\a!) = |0, —1/2). The coupling between these two states 
can be made very large thanks to the Feshbach resonance 
at 834 Gauss. On the other hand, if we choose the final 
state of the Raman process to be \(3) — |1, 1/2) with the 
same nuclear spin component as a, the interaction of \(3) 
with both |a) and \a') is non-resonant, corresponding to 
a low value of the scattering length a a ^p ~ a a >^ ~ 2.5 
nm [21]. This yields a typical scale for the interac- 
tion energy between an atom in (3 and the background 
in a, a' which is smaller than the typical bandwidth, 
and hence negligible. Furthermore, taking typical val- 
ues for the lattice wavelength A ~ 800nm, the atomic 
density p ~ (2/A) 3 and the recoil velocity v ~ h/(MX) 
of atoms in state (3, we evaluate the collision rate to be 
7 C = pav ~ 10 2 s _1 . The Raman detection sequence can 
therefore be performed in a time scale of the order of 
a few milliseconds, yielding an energy resolution in the 
100Hz range. Losses due to inelastic transitions from 
state (3 have a rate ~ 10~ 12 cm 3 s _1 and can be neglected 
on this time scale. 

In summary, we have proposed a Raman spectroscopy 
technique which, analogously to ARPES in solid-state 
physics, is able to probe the one-body Green's function. 
This technique can be used to obtain information on 
the Fermi surface, and on the quasiparticles (or absence 
thereof) of a gas of fermionic atoms, even in strongly- 
correlated states. In the near future, this technique may 
play an important role in the experimental characteriza- 
tion of the novel quantum states of matter that can be 
obtained with ultracold atoms in optical lattices. 
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